An experimental test of volume-equilibration between granular systems 
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Understanding granular and other athermal systems requires the identification of state variables 
which consistently predict their bulk properties. A promising approach has been to draw on the 
techniques of equilibrium statistical mechanics, but to consider alternate conserved quantities in 
place of energy. The Edwards ensemble [l], based on volume conservation, provides a temperature- 
like intensive parameter called compactivity. We present experiments which demonstrate the failure 
of compactivity to equilibrate (via volume-exchange) between a pair of externally-agitated granular 
subsystems with different material properties. Nonetheless, we identify a material-independent 
relationship between the mean and fluctuations of the local packing fraction which forms the basis 
for an equation of state. This relationship defines an intensive parameter that decouples from the 
volume statistics. 
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Granular materials such as sand exhibit a transition 
from solid-like to liquid-like states which cannot be pre- 
dicted from classical thermodynamics due to the separa- 
tion of energy scales: thermal fluctuations are insufficient 
to cause grain rearrangements. However, because the be- 
haviors of granular materials are strongly reminiscent of 
ordinary phases, there have been extensive efforts to for- 
mulate statistical theories as explanations. The Edwards 
proposal [1^ to formulate a volume-based (as opposed to 
energy-based) ensemble of states has provided the basis 
for much work in the field [HHHiailiaElll]. Within 
this approach, a temperature-like parameter called com- 
pactivity plays the role of the central state variable, yet 
fundamental properties such as equilibration have not 
been established. 

Following the Edwards approach, several ensembles 
have been proposed in order to define the associated mi- 
crocanonical entropy. In general, these ensembles are 
grounded in partition functions relying on grain-scale 
spatial tessellations [71 E] which have recently been 
generalized to include stresses [H |TOl [TTJ [121 [E] • each 
case, the temperature- like intensive variable appears in 
the denominator of a Boltzmann-like weight and should 
take the same value in subsystems allowed to exchange 
the conserved quantity. For the Edwards ensemble, the 
microcanonical entropy is S(y) = klogO. where is the 
number of mechanically stable configurations compatible 
with the overall volume V. The associated intensive pa- 
rameter X is known as the compactivity and is defined, 
by analogy with thermodynamics, as ^ = The com- 
pactivity measures how far the system is dilated above 
its most compact state {X = 0). While several mea- 
surements of the compactivity have been performed in 
experiments HI [SJ [6l [9] , none has yet assessed equilibra- 
tion between subsystems. 

In this paper, we present experiments which demon- 
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strate that compactivity does not equilibrate between 
granular subsystems exchanging volume and kept in a 
dynamical stationary state. The two subsystems have 
the same number of grains and differ only in their mate- 
rial properties, but are observed to occupy different total 
volumes. Encouragingly, we find a species-independent, 
one-to-one relationship between the mean local packing 
fraction (j) and its fluctuations. Together, these two re- 
sults provide different values of the compactivity for the 
two subsystems. While the observed constant ratio be- 
tween {(j)) and its fluctuations provides a well-defined in- 
tensive parameter, it cannot provide a prediction for the 
value of the two volumes. This points to the need for an 
ensemble that captures additional aspects of the granu- 
lar material close to its jamming transition, such as the 
coupling between volume and internal forces and/or dy- 
namical features. Such a broader theory must provide 
other intensive quantities which do equilibrate, allowing 
for the establishment of the full equation of state. 

We study the equilibration of two adjacent dense bidis- 
perse layers of disks on an air table, (see Fig. [T]). A pis- 
ton separates the two subsystems, and is constrained to 
move only along the long axis of the air table. On each 
side, the system is prepared with an identical number of 
grains N ^ with the ratio of small to large grains fixed at 
Ns = 2Nl. This ratio suppresses crystallization and pro- 
vides roughly the same volume of large (diameter dL = '^^ 
mm) and small grains {ds = 58 mm). The disks within 
each subsystem have a different set of material proper- 
ties: on Side 1, the particles have restitution coefficient 
ei = 0.51 ± 0.07 and friction coefficient jii = 0.85; on 
Side 2, the particles have 62 = 0.33 ± 0.03 and /i2 = 0-5. 
Both sides utilize standard plastic Petri dishes as the par- 
ticles, with the difference in particle properties achieved 
by encircling the particles on Side 1 with a rubber band. 
Restitution coefficients were measured from isolated bi- 
nary collisions; friction coefficients are nominal values 
from the literature. Because the sides of the particles 
slope inwards, the thickness of the rubber band does not 
significantly change the radius of the particles; the mass 
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FIG. 1: (a) Apparatus schematic (1 m x 2 m), with particles 
drawn to scale. Filled grains are taken from an experimen- 
tal configuration and indicate the size of the imaged region. 
A measures the fractional position of the piston, (b) Sam- 
ple probability density functions of A (and Gaussian fits) for 
N = 183, 192 and 201 (0 = 0.766, 0.789, 0.812). (c) (A) for 
Side 1 as a function of N, with bars representing the standard 
deviation. For some values of the experiment has been 
repeated several times starting from different initial configu- 
rations, in order to evaluate the errors; these measurements 
are shown as smaller symbols. 



of the particles on Side 1 is increased by 7%. 

The aggregate rearranges via an array of sixty electro- 
magnetic bumpers which form the walls of the system. 
These bumpers are triggered pairwise: bumpers facing 
each other in the system fire at the same time in order 
to prevent net momentum and torque injection. Four 
pairs of bumpers are randomly fired every 0.1 second, 
and travel 1 cm into the granular pack; the total time 
during which the bumpers stay in the "active" position 
is ~ 0.095 second. 

To quantify the long-time mobility of the particles, 
we take images at a frequency which is low compared 
to the energy injection timescale, making usual track- 
ing techniques inoperative. Therefore, we have devel- 
oped a tracking method which identifies each particle by 
a unique tag. Each particle is marked with a 3 x 3 array of 
colored dots which encodes two copies of a 4-bit, 4-digit 
identifier, plus an error-correcting bit. The particles are 
located by their circular rims and their identities are es- 
tablished using the tags, allowing adjacent image frames 
to be connected into trajectories. We monitor the posi- 
tions of the piston and the inner 75% of the disks with a 
CCD camera mounted above the apparatus; we obtain a 
minimum of 10^ configurations for each experiment. 

To understand the equilibration of the system on its 
approach to jamming, we perform experiments at in- 
creasing values of while holding all other variables 
constant. Over the studied range N = 183 to 204 




FIG. 2: Diffusion length cr as a function of lag-time r. Thin 
lines are for Side 1; thick lines are for Side 2. Inset: Tra- 
jectories of a single disk (from Side 1) for a duration of 10^ 



(corresponding to ^ = 2 (^j^ + (^2)) ~ 0.768 to 
0.818), the system transitions from a liquid-like to a solid- 
like state. To monitor this transition, we use particle- 
trajectories to co mpute the average diffusi on length, de- 
fined as a (r) = v'(||r;(t + r) -fi{t) where u {t) is 
the vector position of disk i at time t. Fig. [2] shows a plot 
of this quantity as a function of r for three values of 
and for each side. The inset shows typical trajectories at 
corresponding values of A". At A" = 183 (low ^) a parti- 
cle explores a region several particle diameters wide, and 
a saturates at long time scales due to finite system size. 
On the other hand, at A" = 201 (high ^), a exhibits a 
subdiffusive plateau at short r, and caging effects are sig- 
nificant, as can be seen in Fig.|2j This change in mobility 
occurs as (j) approaches and crosses (^rlp- Therefore, as 
the dynamics of the system slows down, we increase the 
duration of the experiment from twenty hours up to fifty- 
five hours for the densest packing in order to sample a 
significant set of configurations. 

The scaled position of the piston, < A < 1, mea- 
sures the macroscopic state of each subsystem. After a 
transient, the probability distribution of A becomes sta- 
tionary and Gaussian, even for the largest A"; sample 
distributions are shown in Fig. flfc. We examine (A) as a 
function of A' and observe that (A) < ^. On average and 
for all values of A", Side 1 occupies less than half of the 
overall volume. We conclude that this systematic devia- 
tion from equal volumes originates from the difference in 
material properties of the disks. Therefore, the density 
of states depends not only on the available space, but 
also on the grains' properties. 

In order to extract a microscopic, "canonical" quantity, 
we analyze the temporal fluctuations of the average local 
(j) over windows of increasing size. For each subsystem, we 
measure (j) over boxes of size L ranging from a few parti- 
cle diameters up to half the system size. We observe that 
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FIG. 3: Variance (^0 ) of the packing fraction measured 
within squares of size L for three values of N] dashed hues 
are ((502)o/L^ 



the variance ((5^^) scales approximately as ^ as can be 
seen on Fig.[3| Such scaling behavior is expected from the 
central limit theorem in the absence of long-ranged spa- 
tial correlations in (j). We can therefore conclude that the 
packing fraction is a self- averaging quantity, with statis- 
tics suitable for a thermodynamic-like analysis, despite 
the relatively small number of particles. We obtain an 
L-independent measure of the variance of (j) by averaging 
= (^02^^2 Q^g^ L > 2ds. 

This analysis provides both the mean packing frac- 
tion {(j)) and its normalized fluctuations {S(j)'^)o defined 

as {5(jP'){L) — ^^^2^. These packing fraction statisitics 
are plotted as a function of N in Fig.|4^,b. The behavior 
of the two sides reproduces our macroscopic measure- 
ment: ((/)) is higher for Side 1. At = 195 the packing 
fractions are equal (((/>) = ^ = 0.798 for both sides), as 
shown in Fig. [4^. This value is in approximate agree- 
ment with independent measurements of random loose 
packing for these grains: (j)^^^ = 0.807 ± 0.010 and 
0^^^ = 0.812 ±0.006. To determine 0^^^ for each of the 
two types of particles, we placed the table at a 0.3° angle, 
and rained down single particles from the high end to the 
low end to create the loosest packing accessible to us. In 
addition, this value of N corresponds to a transition from 
Brownian to caged dynamics (see Fig. [2|. 

Different values of the normalized (p fluctuations {S(j)'^)o 
are observed on the two equilibrated sides, indicating 
that the statistics of this quantity depend on the ma- 
terial properties. Moreover, the fluctuations on each side 
are found to be decreasing functions of A', as shown in 
Fig. |4]3. This can be understood by the fact that as A" 
increases, the amount of free volume to be distributed 
among the grains decreases. Finally, the discrepancy in 
the packing fraction statistics between the two sides gets 
larger on approach to jamming at large 4> (large A"). 

To examine the relevance of these state variables, we 
plot the dependence of {S(j)'^)o on (0) (see Fig.[4|3). Re- 
markably, the data from both sides fall on a master curve: 
the fluctuations in are uniquely determined by the local 
free volume, and are insensitive to the material proper- 
ties of the grains. Within the explored range of A^, this 
master curve is approximately linear and the extrapola- 
tion would intersect (^(/>^)o = at {(j)) = 0.845, which is 




FIG. 4: (a) Packing fraction (0) and (b) normalized packing 
fraction fluctuations ((50^ )o as a function of number of parti- 
cles on each side, (c) Data from parts (a) and (b) combined 
on a single plot without regard for N. Ellipse encloses two 
points obtained from Side 1 (red) and Side 2 (black) for a 
single run with N = 186. Dashed line is a linear fit to data 
from both sides, showing an intercept at = 0.848. 



compatible with the rigidity transition reported in other 
bidimensional systems [Ml |T5l [16] . 

In order to measure the compactivity, we follow the 
method of Nowak et al. [2^ We will consider that each 
side is a compactivity reference, or compactostat^ and use 
the microscopic information as a compact ometer. Within 
this canonical framework, we start from the Edwards par- 
tition function Z = e~^^^^/^ where the summation 
runs over the set C of mechanically stable configurations 
of the assembly, and X is the compactivity. In our sys- 
tem, configurations above (j)^^^ are mechanically stable 
and thus the ensemble explored by the rearrangements 
closely matches the Edwards ensemble. By analogy with 
the standard canonical ensemble, successive derivatives 
of In Z provide higher moments of the volume distribu- 
tion. In particular, the mean volume and its fluctuations 
are related by 



(1) 



We integrate Eq. [T] from a reference state with known 
values (Xref,Kef) to the current state (X, F), yielding 
the relation 



dV 



(2) 



In our experiments, the ^-statistics are observed to be 
identical for the two subsystems at A" = 195. Therefore, 
we can use this as a reference state with Xref and V^^i 
equal for both sides. In Fig.jiJ:, we observe that the inte- 
grand in Eq. [2]is the same for both sides. Therefore, we 
can perform the integral from A'ref = 195 to any value 
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of and obtain the compactivity measured relative to 
Xref . We are particularly interested in the jammed states 
{N > 195): performing the integral over the same inte- 
grand to two different values of V yields two different 
values of X. Therefore, Xi ^ X2 even though the sub- 
systems are equilibrated (in steady state). 

This conclusion relies on the unexpected one-to-one re- 
lationship between the average packing fraction and its 
local fluctuations. This relationship should be predicted 
from an as-yet unknown equation of state. It is tempting 

to define a variable Y = ^^^y^ which is intensive, (j)- 
independent, and takes the same value in both (all) sub- 
systems. However, the value of F, which we conjecture 
to be related to the energy injection and/or dissipation 
rates in the system, is not enough to specify the equili- 
brated volumes. This suggests that another parameter 
is needed to fully characterize the state of a dynamically 
evolving dense granular pack, much as both temperature 
and pressure equilibration are needed to solve the equiv- 
alent problem within classical thermodynamics. 

It is instructive to speculate about the nature of the 
missing information. First, the Edwards ensemble relies 
on an assumption that configurations occupying the same 
volume are equiprobable. Recent simulations [17 on fric- 
tionless particles have challenged the validity of this as- 
sumption. However, even for unequal probabilities, the 
integral in Eq. [2] remains valid [6 , and thus the mea- 
sured compactivity would still fail to equilibrate. Second, 
higher ja or lower e might provide for more efficient energy 
dissipation in the bulk and couple the grains differently 
to the energy injection at the boundaries. Measuring 
energy transfer rates might prove useful to decipher the 
statistics of such a system and studies of grain ki- 
netics would allow further investigation of such a mecha- 
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nism. An approach based on a fluctuation/response ratio 
[19 provides an additional notion of effective tempera- 
ture closely related to these dynamical features. Finally, 
and importantly, intermittent force chains are particu- 
larly abundant in higher-density states [20 and likely 
couple the microscopic stress state to the macroscopic 
volume. Recent proposals [lOj HTJ [T2j |13] to include con- 
tact force or stress statistics in defining intensive statis- 
tical parameters might provide the solution. 

These experiments have examined microscopic and 
macroscopic jammed matter statistics in light of existing 
theories. Remarkably, a grain-independent relationship 
was found between packing fraction average and fluctu- 
ations which provides a well-defined state variable, al- 
though insufficient to fully characterize the state of the 
dynamic packings; further investigations are needed in 
order to establish the generality of this observation. Our 
measurements also highlight the important role played 
by the static cj)^^^ in the dynamics of a granular sys- 
tem. The connection between the jamming point 
(/)^, and granular glassy dynamics [21] should be further 
investigated. Finally, we observed that subsystems with 
different material properties do not in general share vol- 
ume equally. This underlies the crucial role played ma- 
terial properties such as friction in determining the state 
of a granular packing. 
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